import qgis
from qgis.analysis import QgsRasterCalculatorEntry, QgsRasterCalculator
from qgis.core import QgsProject
from qgis.analysis import QgsZonalStatistics
from PyQt4.QtCore import *
from osgeo import gdal, ogr,osr
import sys
import processing
from os.path import join

def extraer_media(vector,p_raster,campo):
    zoneStat= QgsZonalStatistics(vector,p_raster,campo,1,QgsZonalStatistics.Mean)
    zoneStat.calculateStatistics(None)

def copia(vector,salida):
    crs=vector.crs()
    clonarv = QgsVectorFileWriter.writeAsVectorFormat(vector,salida,'utf-8',crs,"ESRI Shapefile")
    

sig_path= "C:/Dropbox (LANCIS)/SIG/desarrollo/sig_megadapt/procesamiento/"
p_ageb_o=sig_path+"encharcamientos_escorrentias/precipitacion_15_16/agebs.shp"
p_ageb_prec_16=sig_path+"encharcamientos_escorrentias/precipitacion_15_16/salida/tp_ageb_prec_16_v1.shp"
v_ageb=QgsVectorLayer(p_ageb_o,"ageb","ogr")
copia(v_ageb,p_ageb_prec_16)
v_ageb_prec_15=QgsVectorLayer(p_ageb_prec_16,"ageb_prec","ogr")
p_prec=sig_path+"encharcamientos_escorrentias/precipitacion_15_16/salida/prec_16.tif"
campo ="prec"

extraer_media(v_ageb_prec_15,p_prec,campo)
print "proceso terminado"

import qgis
from qgis.analysis import *
from qgis.core import QgsProject
from PyQt4.QtCore import *
from osgeo import gdal, ogr,osr
import sys
import processing
from os.path import join


sig_path= "C:/Dropbox (LANCIS)/SIG/desarrollo/sig_megadapt/procesamiento/"
p_ageb_prec_16=sig_path+"encharcamientos_escorrentias/precipitacion_15_16/salida/tp_ageb_prec_15_v1.shp"
vector=QgsVectorLayer(p_ageb_prec_16,"ageb_prec","ogr")
p_salida=sig_path+"encharcamientos_escorrentias/precipitacion_15_16/salida/tp_ageb_prec_15_v2.shp"
campo="precmean"
n_c_norm="prec_n"

def normalizar(vector,p_salida,campo,n_c_norm):
    #crea una copia 
    crs=vector.crs()
    clonarv = QgsVectorFileWriter.writeAsVectorFormat(vector,p_salida,'utf-8',crs,"ESRI Shapefile")
    
    lista=[]
    vector_s=QgsVectorLayer(p_salida,"salida","ogr")
    for f in vector_s.getFeatures():
        lista.append(f[campo])
        
    minimo = min(lista)
    maximo = max(lista)
  
    vector_s.dataProvider().addAttributes([QgsField(n_c_norm,QVariant.Double)])
    print maximo,minimo 
    vector_s.startEditing()
    
    for f in vector_s.getFeatures():
        normalizar=(f[campo]-minimo)/(maximo - minimo)
        f[n_c_norm]=normalizar
        
        vector_s.updateFeature(f)
    
    vector_s.commitChanges()
    


normalizar(vector,p_salida,campo,n_c_norm)

print ("proceso terminado")